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Abstract 

Structural disorder is present in almost all experimental measurements of electronic transport through sin- 
gle molecules or molecular wires. To assess its influence on the conductance is computationally demanding, 
because a large number of conformations must be considered. Here we analyze an approximate recursive 
layer Green function approach for the ballistic transport through quasi one-dimensional nano-junctions. We 
find a rapid convergence of the method with its control parameter, the layer thickness, and good agreement 
with existing experimental and theoretical data. Because the computational effort rises only linearly with 
system size, this method permits treatment of very large systems. We investigate the conductance of gold- 
and silver wires of different sizes and conformations. For weak electrode disorder and imperfect coupling 
between electrode and wire we find conductance variations of approximately 20%. Overall we find the 
conductance of silver junctions well described by the immediate vicinity of narrowest point in the junction, 
a result that may explain the observation of well-conserved conductance plateaus in recent experiments on 
silver junctions. In an application to flexible oligophene wires, we find that strongly distorted conformations 
that are sterically forbidden at zero temperature, contribute significantly to the observed average zero-bias 
conductance of the molecular wire. 

PACS numbers: 



1 



I. INTRODUCTION 



The study of electron transport properties at the molecular or even atomic scale has generated 

nnnn 

many striking insights in the last decade 111,121,131,141]. Many features of the measured transport char- 



acteristics of molecular wires in breakj unctions or AFM/STM setups [15|,16D could be explained on 
the basis of the Landauer formalism. There has been tremendous progress in the characterization of 



metallic nano- wires O, l8|, l9|, 



3, 



llh . organic nano- wires u2. 



ISh and nanotubes 



16D. Un- 



fortunately, due to the complexity of the experimental setup, most of these measurements have no 
control over the details of the electrode geometry, which results in a spectrum of IV-characteristics. 
In the face of these and other uncertainties theoretical models can help understand and explain the 

n 

experimental data [1170 ). These calculations are challenging because the electronic structure of the 
central cluster of the junction („extended molecule"), comprising the molecular wire and some 
part of the electrode fragments, must be computed in atomistic models. Because of the high cost 
of accurate electronic structure calculations, IV-characteristics are often computed for idealized 
geometries and electrode configurations lllSn . 

Due to computational limitations most theoretical studies of molecular junctions have em- 
ployed idealized junction geometries, including placement and contact of the molecule with re- 
spect to the electrode, as well as the electrode geometry. Many present day fabrication techniques 
for nanojunctions entail strong deformations of the system, e.g. rupture of electrode wires (in 
break-junctions) and subsequent contacting procedures, which make it unlikely that ideal junction 
geometries are ever realized in practice. Despite the use of accurate electronic structure methods, 
only semi-quantitative agreement with the measured conductance has been obtained in a number 
of studies to date. Imperfect electrode geometries and contacts are likely to contribute to this 

n 

widely discussed discrepancy HlOn . but are difficult to assess experimentally. 

The effects of structural disorder, of impurities and non-ideal electrode geometries are diffi- 
cult to asses with transport methods that relies on computationally demanding ab-initio electronic 

structure theory, which often scale as 0{N^) with the system size. One possibility to take such 

I I 

effects into account is to develop efficient approximations [20], which allow for inclusion of larger 
electrode fragments in electronic structure calculation of the „extended molecule" region. While 
such approximate methods may not predict the absolute value of the conductance quantitatively, 
they may help explain differential features that may originate from the existence of a large struc- 
tural ensemble in the junction geometry. Such methods are also useful to investigate the coherent 
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transport properties in very large, complex molecules, for which high-accuracy calculations re- 
main infeasible to date. 

I — I 

In previous work tight binding models were successfully applied by Pauly et ah |!21] to inves- 
tigate conductance histograms and strain forces of Ag, Pt and ferromagnetic Ni break junctions. 
Additionally Hafner et al. [l22n studied the influence of the absence of magnetic domains in the 
ferromagnetic 3 J materials Fe, Co and Ni. In the case of metal-molecule-metal junctions based on 
oligophenylene wires the length-dependent conductance and thermopower as well as elastic and 
photo-assisted transport properties, a ;r-orbital tight binding model has given intriguing insights 



In this paper we investigate convergence of a recursive layer Green's function (RLGF) 
approach bsn that scales linearly with the system size. This approach permits the treatment of 



very large systems comprising more than a thousand atoms, e.g. DNA ll26h or the atomic transis- 
tor as well as computations of ensembles of thousands of structural models for disordered 
systems. We first analyze the truncation effects on the accuracy of conductance calculations for 



and 



a several systems, where much experimental and theoretical data is available. Gold [|28|, 
silver |[nl.]2lll molecular wires are among the most widely studied systems in molecular elec- 
tronics. We find that even for metallic systems, where the electronic wave functions are most 
extended, the conductance converges rapidly with the truncation cutoff parameter and converges 



2311 . at a fraction of 



to values in good agreement with experiment and accurate prior studies II18L 
their computational cost. 

Next we apply this method to novel, computationally challenging applications on transport 
through molecular wires. We investigated different scenarios for static and dynamic disorder in 



metal[|7La,l9L 



y 



12 



21 



28,3] 



28LI22D and metal-organic molecular lll3[] wires. In an application on 



the influence of a non- vanishing tilt- and twinning-angle between the source/drain electrode tips 
and find a variation of up to 20% in the zero bias conductance. Additionally we probe the effect 
of randomly introduced surface vacancies on the electrode tips, which can lower the conductance 
by up to 30%. 

We then turn to silver nano-clusters, which are promising systems for applications as electronic 
materials and for surface nano-structuring. While their structural properties have been extensively 



investigated 



3li 



32L 



330, much less is known about their transport properties. We present 



the first theoretical study of electronic transport through large clusters (260 atoms) and study the 
dependence of the transport properties on the structure and size of the nano-cluster. 
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While these examples concern mononuclear systems, the presence of metal-organic junctions 
significantly complicates the physics of nanoscale transport. Charge transfer, differences in the ac- 
curacy of the description of electronic structure of the electrode and the organic molecule, as well 



35h . We have inves- 



as other factors, complicate the theoretical description of such systems [134|, 
tigated the effect of dynamic disorder on the zero-bias conductance of oligophenylene molecular 
wires as a function of their length and geometry. We compare the total transmission function 
obtained by the layer approximation with the full device calculation and find a good quantitative 
agreement at a fraction of the computational cost. 

Recent investigations have focused on the importance of even small scale geometric fluctua- 
tions on molecular transport. Here we extend this work to the first investigation of large-scale 
conformational change. We generate a large thermodynamically relevant ensemble of, which in- 
clude complete ring-flips of the molecules, comprising thousands of conformations, to discuss the 
dependence of the conductance change on the torsion angle fluctuations. We find this system to be 
a striking example, where large scale structural fluctuations significantly affect the conductance. 
Configurations, which are sterically forbidden at low temperature, contribute significantly to the 
overall time-average of the zero-bias conductance. 

II. METHODS 

The electronic conductance of nanoscale junctions is often analyzed in the framework of Lan- 



dauer theory I37n . which divides the system conceptually into idealized electrodes and a central 
region called „extended molecule", that contains the nanojunction of interest (see Figure [I^). Us- 
ing Landauer's formula, the current is computed as: 

m = dEx{E,V)[fL{E)-fR{E)] (1) 

where T denotes the total transmission function. In order to evaluate the transmission function of 
the extended molecule, we must compute its electronic structure in an adequate atomistic model 



I 38h . Equation (1) stipulates the existence of perfect Fermi-seas at the appropriate chemical poten- 
tial of the electrodes. For these quantities to have meaning, electrode fragments must be included 
into of sufficient size into the extended molecule region, so that charge transfer and image-charge 



effects can be taken into account 



350 . The extended molecule for which the transmission 



must be computed may thus be much larger than the physical object which is placed between the 
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electrodes. In particular for homonuclear systems, such as metallic wires, there is no clear-cut 
distinction between electrode and "system" at all. To compute the current, we consider the total 
charge of the system on the right side of an arbitrary division of the system into a right (R) and left 
(L) part, respectively (see Figure 1(b)): 

P = -^ L (2) 

/grhs 

where cj creates an electron in orbital / and the sum runs over all orbitals on the right hand side. 
The current / can then be obtained in perturbation theory from the equation of motion ihl — 
[//,p] = —iehvz, where 

i^^z = L L ~ (3) 

ieL jeR 



defines the velocity operator Il25h . In the expression of the electron velocity operator the indices 
/ and j run over the orbitals in the left (L) and right (R) parts of the system, respectively, is 
proportional to the transfer matrix element tij times the electron propagation, expressed by the 
particle creation and annihilation operators cj and c/, respectively. Calculation of the charge flux 
at any arbitrarily chosen interface is thus sufficient to obtain the total current / = —e{vz) through 
the system. 



380 we can evaluate the coherent zero bias transmission T at T = 



Using Kubo's formula 
as 

T = Tr [(//iv^)Im W{E) {mv^)lm W{E)] . (4) 

The imaginary part of the Green's function Im ^{E) is related to the advanced/retarded Green's 
functions = [^^^ by 

Im W{E) = ^ [W'{E) - W%E)] . (5) 

The computation of in general requires the inversion of the Hamilton operator of the extended 
molecule of the system. The computational effort of this calculation grows rapidly with the system 
size (i.e. the total number of orbitals), restricting the size of the system and the complexity of the 
electrode fragments that can be included in the extended molecule. 

For this reason, we divide the extended molecule into a set of principal layers perpendicular 
to the current flow direction and describe the electrons contained in one principle layer with a 
block Hamiltonian matrix. This introduces a truncation parameter which was found to be well 
converging with an increasing layer thickness (see results section A). We take only these overlap 
matrix elements into account, which belong to nearest neighboring layers and obtain a Hamiltonian 



matrix of band-diagonal type. The main advantage of this approach is that the computational effort 
of the matrix inversion scales linearly with the system size. 

We employ this feature in the evaluation of the transmission function using a Recursive Green's 
function scheme ll25n . This method is based on an iteration over the principal layers from the 
left end to the right end of the extended molecule. For every layer we calculate the Green's 
function by inverting the corresponding block Hamiltonian and take the influence of the semi 
infinite right region via self energies into account. The iteration ends up with the calculation of 
the Green's function of the left most principal layer, which is plugged together with the surface 
Green's function and the corresponding coupling matrixes into Kubos formula for conductance. 

The whole method is based on the assumption of non-interacting electrons permitting the use of 
a tight binding like model Hamiltonian (extended Hiickel), that allows for an efficient computation 
of the electronic structure during a dynamic process. More details of the numerical methods are 
described in the appendix. 

III. RESULTS 

A. Convergence of the layer approximation 

Metallic nanowires have been among the first and most widely studied systems in molecu- 



lar electronics uJi 



211, 



290. Metallic systems are often most challenging for linear- scaling 
electronic structure methods, because the electronic wave functions are extended. To provide a 
stringent test for our „local" approximation, we have investigated the convergence of the layer 
approximation for two representative examples, namely gold- and silver- wires respectively. 

We begin the investigation by dividing silver and gold model junctions into a set of „principal 
layers" with increasing thickness w. In order to establish the convergence of the method for large 
systems, the test geometry has to be of sufficient length. Here we investigate junctions of 45.2 A 
length in z-direction, containing 388 silver or gold atoms with a nearest neighbor distance of 2.88 A 
in both metals Q,!^]- We prepare the electrodes as perfect fcc-clusters, which narrow towards 
the center to form a single-atom point contact at their tips, generating a dimer structure which 
permits a current flow in the crystallographic [111] direction. The extended molecule region and 
the layer divisions are illustrated in Fig. [2^. The bulk electrodes are designated by the two larger 
layers on each side of the system. 
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We calculate the conductance for varying widths w = 1, . . . , 18 d[m] of the „principal layers" 
(Fig. [2^), i.e. the full length of the extended molecule region. Figure [2b shows that the con- 
ductance of the silver and gold model junction as a function of the principal layer thickness w 
converges rapidly to the experimental value. The same holds true for the junction conformations 
labeled Ag 1 - Ag 4 and Au 1 - Au 4, which have a minimum cross section of 1-4 atoms, respec- 
tively. Furthermore, we calculated the convergence of the conductance for longer wire geometries 
(Fig. [it), constructed by sequentially introducing additional atoms into the point contact at the 
minimum cross section. Thereby we obtain silver and gold junctions of 20, 24, 28 and 32 atomic 
layers in the z-direction. 

For w = 1 the conductance is significantly underestimated to approximately 0.5 Gq for both 
metals, indicating that hopping processes across distances larger than the interatomic distance are 
important. For all choices of the layer thickness with w > 1 the conductance has converged to 
nearly the experimental values. For w = 3 the layer division retains the symmetry of the [111] 
crystal stacking order "ABCABC..." in fcc-latices. We investigate the convergence in more detail 
at the level of the transmission in Fig[3l which shows the total transmission function t{E) of the 
geometry shown in Fig. [2^ over an energy interval [Ef — 6 eV, E"/? + 6 eV] . Again we find that all 
curves for w > 2 agree well with one-another. 

In order to demonstrate the efficiency of this method we compare the computation time of the 
transmission curves shown in Fig 2. With a resolution of AE" = 10 meV the transmission of the 
system divided into 1, 18 layers required 1368, 594, 429, 336, and 294 seconds, respectively. 
Using this approximation, e. g. with 6 principal layers, that takes 31% of the time of the „full- 
device" calculation, while increasing only a neglectable los of accuracy. 

B. Imperfect Electrode Geometries 

There are two obvious parameters which define the junction geometry with respect to the elec- 
trodes that are presently not under experimental control: the tilt and twinning angles of the two 
electrode fragments with respect to one another. In order to investigate the dependence of the con- 
ductance on these parameters, we have prepared a ideal fee- silver junction with 224 atoms as in 
the previous section and varied the tilt (a = 0, ...,70 deg) and twinning (jS = 0, ...,60 deg) angle, 
as shown in the insets of Fig. |4l Increasing the tilt-angle a from to 20 degrees leads to a slight 
increase of the conductance by circa 0.05 Gq which can be explained by the influence of interfer- 
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ence effects, which strongly depend on small changes of the atomic positions. A further increase 
of a from 20 to 70 degrees results in a decreasing conductance by 0.15 Gq, which corresponds to 
the loss of crystal symmetry across the junction. In contrast, twinning the electrode from to 60 
degrees leads only to a minor change in the conductance of 4 • 10~^ Gq. 

When an electrode is manufactured in a break junction or generated by contacting the tip of an 
AFM/STM, it is very unlikely that the perfect lattice geometries with perfect surfaces along the 
crystalline axis, assumed in nearly all theoretical investigations, are realized in practice. According 
to all models of electronic transport, each surface defect creates an additional scattering center 
that may impede coherent transport through the junction. Imperfections in the geometry of the 
electrode tips will thus influence the ballistic transport. On the other hand, we have seen in the 



previous section, in agreement with many prior studies 11 17[ 



2ll, 



28,29|] 



28,29n, that the conductance of the 



junction is mostly determined by its most narrow region 114 Ifl . 

In order to estimate the significance of tip disorder we have therefore prepared a perfect junc- 
tion as above and then randomly removed atoms from the surface of the electrode in the vicinity 
of the contact point. The number of silver atoms in the extended molecule region is systematically 
decreased by removing 28 atoms at randomly chosen surface positions. To maintain coherent 
transport, the two central atoms were never removed. Every junction geometry with n • 28 vacan- 
cies = 1, ...,5) was generated 500 times, with randomly chosen vacancy positions. For each 
conformation we computed and subsequently averaged the conductance. The calculations were 
performed using a principal layer thickness parameter w = 3. 

Samples of the resulting junction conformations are shown in Fig. [5] with the corresponding 
averaged conductance value in units of Gq, respectively. In addition, the total number of surface 
vacancies on the current junction geometry is given below the conductance values. Figure [6] shows 
the resulting conductance values averaged over conformations with equal number of surface im- 
purities. Creating 140 vacancies, which is half of the initial number of atoms, reduces the total 
conductance by 40% in average. The rapidly increasing size of the error-bars indicates that the 
change in the conductance depends strongly on their positions. 



C. Silver Nanoclusters 



Recent experiments of silver junctions yjllllt] suggest a strong stability of the observed zero- 
bias conductance in electrochemically grown silver junctions. While this effect was locally 
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explained ll27h by the selection of specific contact geometries, the overall shape of the silver con- 
tacts is likely to vary strongly from one realization of the next. To assess the effect of these 
large-scale geometric differences, we have prepared locally similar, but globally different junction 
geometries by placing silver clusters of various size in different orientations on a perfect surface 
and then contacting the tip of the cluster with an "ideal" junction. For this purpose we use the 



optimized cluster geometries from Ref. 13 in . which where generated by Monte Carlo mini- 



mization and the modified dynamic lattice search method. Figure |7^ shows the top view of the 
studied silver clusters with 5, 7, 180, 220, and 260 atoms and decahedron (m-Dh) core symmetry. 
As illustrated in Fig. Uh we consider the metal clusters attached to a silver substrate layer of the 
crystallographic [111] direction. The second electrode is realized by a pyramidal tip on top of the 
nano-cluster similar to a STM-setup. 

We optimized the position of the silver nano-cluster on the substrate using a Metropolis Monte- 
Carlo technique combined with the semi-empirical Gupta potential for the silver atoms as de- 
scribed above. During the simulation the silver cluster is treated as a rigid body, so only trans- 
lations and rotations of the cluster are allowed - structural rearrangements insight the cluster are 
forbidden. The metal cluster surface consists of a set of [111], [110] and [100] facets. The mini- 
mum of the potential energy is reached, if the system is arranged such that the largest [111] facet 
(which is always the largest subsurface in the present cases) and the [1 1 1] substrate layer are fac- 
ing each other. The top electrode is assumed to point directly on one arbitrarily chosen silver atom 
on the cluster surface. 

Figure!?]: shows the total transmission function of the clusters Ag5, ... , Ag26o- The conductance 
of the systems is given by the average of the transmission over a small interval around the Fermi 
energy [Ef — A,£'/ + A] with A = 50 meV. For the clusters with 5, 7, 180, 220, and 260 atoms we 
find conductance values of 1.10, 1.08, 1.15, 1.17, and 1.17 Gq, respectively, which means, that the 
conductance is less effected by the size of the nano-cluster and depends more on the point contact 
to the second electrode. This observation may explain the observed stability of the experiment: 
While reconstruction of the junction geometry assures the selection of a specific local geometry, 
the overall conductance depends only very little on the global shape of the clusters forming the 
contact. This result is also in good agreement with an earlier study using a simplified model 
based on random matrix theory. 

We also note that an irregular fluctuation of the transmission as a function of energy is ob- 
served, which increases with the cluster size. Such fluctuations can be conceptually explained by 
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the interference of the incident electron waves with waves scattered repeatedly in the extended 
molecule region containing the sliver nano-cluster and the electrode tip. Increasing the cluster size 
permits scattering processes of increasing order which leads to conductance fluctuations of higher 
frequency. An analysis of the average energy spacing of the extrema of the transmission (which 
may be measured by applying a gate voltage) can help to estimate the size of the backscattering 
region. 



D. Thermal fluctuations of the conductance of oligophenylene wires 



So far we have discussed only mononuclear extended molecule regions. It is well known that 
the presence of metal-organic interfaces complicates electronic structure and, as a result, electronic 
transport calculations. We have therefore studied the coherent conductance of phenyle-di-thiol 
(PDT), a „Drosophila" of molecular electronics. Because this molecule was studied extensively 



in the past 



JJrosopmi 



42|,|43, 



45 



460, it allows for a comparison of the RLGF approach with 



experiment and various other levels of theory. We investigate the transmission of oligophenylene 
molecules of varying lengths, which lend themselves nicely to an investigation of the layer ap- 
proximation in an organic, semi-conducting system. The structure of the molecules suggests a 
natural introduction of J_ayers in terms of single phenyl-ring units, similar to the layers introduced 
by base-pairs in DNA MTfl . 

Figure [8^ illustrates the oligiphenylene molecules covalently bound to Aui9-clusters, using the 
same notation as for the conformations in Ref. kSfl . The molecule is connected to the Au elec- 
trodes at both sides by a symmetric covalent bond of a sulfur atom to three Au atoms. In the liter- 
ature this bonding situation is referred to as the hollow position ll48n . The electrode clusters where 
constructed from fee lattices as above, while the geometry of the phenyl wires was optimized using 



500 . As can be 



density functional theory (DFT) in the local density approximation (LDA) [|48|, 
seen in Fig. [8^ there is a non- vanishing tilt-angle between the phenylene-rings, which varies be- 
tween 33.7 and 34.5 degrees due to the interplay of steric repulsion and ;r-conjugation of adjacent 
rings. A detailed investigation of the influence of (conjugation induced) tilting on the coherent 
transport properties in biphenyl-derived dithiols was recently given in Ref. Isij . 

In Fig. [8h the total transmission is shown as a function of energy for the molecular junctions 
above, once with (dashed line) and once without the layer approximation (solid line). Due to 
the neglect of several overlap matrix elements the transmission and conductance obtained with the 
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layer approximation is below the full device transmission (see tableU). The measured conductances 
of oligophenylene wires with amine end groups indicate even lower conductance values (also given 
in table |I]) which may arise from the differences in the coupling to the electrodes. With increasing 
length of the phenyl wire the transmission gap decreases from 3.89 eV to 2.56 eV. The equidistant 
transmission at the Fermi energy of the different molecular wires indicates the correct exponential 

nnnn 

decrease of the conductance with linear increasing wire length [124J,|52|,|53L154J]. The proportionality 
of the conductance decay G/Gq ^ e~l^^ of the oligophenylene wires once with and once without 
the layer approximation is shown as inset in Fig. [8b. In both cases we obtain a decay coefficient 
jS = 1.47 which is close to the experimental value ^exp. = 1.5 reported in [l54ll for amine end 
groups. Nevertheless, the comparison between theory and experiment remains difficult because of 
the different end groups used. 

Next we investigate the influence of thermally induced molecular vibrations on the coherent 
transport properties of an Au-h-R4 wire. For the simulation of the dynamics of the system we use 
of the AMBER 8 molecular dynamics package 115 511 . which employs the well established GAFF- 
forcefield and a Langevin thermostat method to model temperature. Assuming fixed gold atoms 
of the electrodes we simulate the evolution of the system at 300 K for 10 ps in time steps of 2 
fs. Every 10-th time step a snapshot of the conformation is taken as input for the conductance 
calculation generating 500 conformations for analysis. For each conformation we calculate the 
zero-bias conductance. Within the simulation period we find repeated conductance fluctuations by 
more than an order of magnitude. 

Recent investigations have already focused on the influence of intramolecular vibrations on the 
conductance ll27n . Here, we find an interesting model system where thermal fluctuations induce 
large-scale conformational change. The conductance of a conformation correlates highly with 
its "planarity", because the fully planar conformation leads to a strong overlap of the ;r-orbitals, 
which in turn increases the transmission. However, such planar configurations are forbidden at 
zero temperature because of steric repulsion of the hydrogen atoms emanating from the rings. 
Figure [9^ shows the fluctuation of the torsion angles , (j)2 and (j)^ between the ring-units occurring 
in Au-h-R4, respectively. All three torsion angles fluctuate strongly around an equilibrium average 
of 33.9 degrees, which agrees well with the experimentally observed equilibrium value of 34 
degrees. The histogram Fig. [9b shows that the outer torsion angles ^\ and (/)3 have slightly broader 
distributions than which may be caused by the lower potential energy barrier at the electrodes. 
The average conductance over 10 ps arises as an average of strongly fluctuating instantaneous 
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values, as illustrated in Figure [9]:. 

In the course of the the 10 ps simulation we find 6 "near-planar" configurations of adjacent rings 
and 2 "near-planar" conformations of all three rings. To characterize this geometric feature we 
compute the average of the absolute torsion angles = | + |(/)2| + |(/)3|), which has 4 minima 
(with < 20 deg) and 3 maxima (with > 45 deg) in the observation period. As illustrated in 
Fig [It the conductance has associated minima where ^ is maximal, e. g. at 2.2, 4.3, and 9.4 
ps. Correspondingly the highest conductance values are obtained, if ^ is small, e. g. at 4.6, 5.5, 
and 9.8 ps. This effect is also illustrated in Fig.[9]i, which gives a a higher time-resolution of the 
grey regions of Fig. [9^ and [9t. This analysis shows that the experimentally relevant conductance 
at room temperature arises as an average of thermally excited high-conductance conformations, 
which are forbidden at zero temperature. 



IV. CONCLUSIONS AND DISCUSSION 



Most experimental realizations of nanoscale junctions will contain some degree of structural 
disorder, which is difficult to assess in situ experimentally. The importance of thermal fluctuations 



on the local geometry was recent 
with experimental observation!^ 



y modeled for benzene molecular junctions in good agreement 



20h . Here we use a recursive layer Green's function approach 



for the ballistic electronic transport through several disordered nanoj unctions. 

We have investigated the influence of imperfect contact geometries on the conductance by 
studying many different possible realizations of silver and gold nanowires. Tilting the junction 
electrodes up to 60 degrees reduces the conductance by 20%, while twinning the electrodes leaves 
the conductance nearly unchanged. We found that the introduction of up to 50% surface vacancies 
leads to only small variates of the ballistic transport properties of silver contacts, as long as the 
immediate vicinity narrowest point of the junction was not affected. This analysis was supported 
by a study of the effects of global conformational change in silver junctions for locally conserved 
junction geometries. Here we find that variations in the global cluster geometry ranging from 5- 
260 atoms have only a weak effect on the zero-bias conduction of junctions with locally conserved 
geometry. 

Extending earlier work on local fluctuations we have investigated the impact of thermally in- 
duced large-scale geometric change on the conductance of oligophenylene wires. Averaging the 
conductance over 500 conformations obtained from a 10 ps molecular dynamics simulation at 
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room temperature, we find temporal conductance fluctuations of more than an order of magnitude. 
The average conductance, which agrees well with the experimentally observed value[|52[], results 
from high-conductivity conformations that are sterically unfavorable at zero temperature. The un- 
derlying advances in the modeling approach thus permit the detailed characterization of disorder 
effects, which present in almost all experimental realizations, on the conductance of molecular 
wires. 
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V. APPENDIX 

A. Layer approximation 

The key approximation of the recursive-layer Green's function approach is the division of the 
quasi one-dimensional device region into „principal layers" (see Fig. 1(b)) perpendicular to the 
current flow bSfl . In the layer- approximation all Hamiltonian overlap matrix elements hij between 
atoms which are separated by more than one principal layer are set to zero. Therefore, the full 
Hamiltonian H of the device region has a block-diagonal form, with = {hij} representing the 
Hamiltonian of the „principal layer" k. 

\ 



H 



Hi 


H\2 







Hi 


H23 





H32 


H3 











(6) 
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This approximation seems natural for a number of systems, e.g. polymers such as DNA, where 
the matrix elements between neighboring base-pairs decay rapidly with distance. Overlap matrix 
elements, as well as Hamiltonian matrix elements of localized atom-centered basis sets decay very 
quickly with the nuclear distance. Once such a set of principal layers is chosen, we can calculate 
the zero-bias conductance applying a recursive Green's function method. In the limit of small bias 
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voltages Verges ll25h showed for a rectangular lattice of sites, that the conductance can be evalu- 
ated by a recursive calculation of the „column" Green's function, starting at the right lead iterating 
to the leftmost site column of the device region. The knowledge of the charge flux in the contact 
between the left most column and the left lead is sufficient to evaluate the total current of the sys- 
tem (due to flux conservation). Here we generalize this scheme to a non-rectangular, all- valence 
electron model to permit efficient material specific conductance calculations in arbitrary junctions. 
Since molecular wires are "quasi one-dimensional" by construction, these systems are always nat- 
urally dividable into principal layers in the transport direction. The layer width, which interpolates 
between the full system and single-atom layers, is a free parameter of the approximation, which 
must be chosen to reflect the properties of the system. 



B. Recursive layer Green's function method 

Let us in the following denote as layer and N+ltht entire left and right electrode fragment, 
respectively (see Fig. 1(b)). The computation of the Green's function starts at the rightmost layer 
of the central region of the system, containing the rightmost device layer N and the semi-infinite 
right lead layer N+ 1. Its electronic structure is reflected in the retarded Green's function matrix 

^ I ESn-Hn ESn,n+\-Hn^n+\\ 

\ESn^ 1 ,7V — 1 ESn^ 1 — Hj^^ 1 J 

which is more conveniently expressed in terms of the layer- self-energies Z^. The self-energy Ejy 



of the right contact can be computed from the retarded surface Green's function g^{E) jlSu as 

= {ESn^i^n — Hn+i,n){ESn^i—H]\j^i) ^ 

X {ESn^n^i —Hn^n^i) (8) 
^ Ta.+ 1,A.(^) gjv+l(^) 4+i,A^(^) (9) 

Here T/j denotes the coupling matrix of the layers / and j. We can calculate the Green's function 
and self-energy of every principal layer k recursively, using the relations: 



gl{E) = (ESk-Hk-mE))-' (10) 



i:l_,{E) = Tk-i,k{E) gi(E) (11) 



with k = N 2. 
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In a layered system, we have chosen the interface such that the velocity operator has non- vanishing 
terms only for orbitals connecting the left-most electrode layer (layer 0) with orbitals in layer 1 . 
To compute the GF occurring in equation ([5]) we thus need only the retarded Green's function of 
the system comprising layer and layer 1 

^\E) = -//oi -^(^)]-\ (12) 

which is easily computed from the right- and left-lead self energies Z-* (/ = 0, 1). 

The final ingredient to the calculation is the computation of material- specific electrode Green's 
functions for the left and the right reservoir. In order to avoid problems that may arise from an 
incommensurate decription of the electrodes and the extended molecule, we compute the electrode 
Green's functions using the same parametrization of the Hamiltonian. We assume a semi-infinite 
fcc-lattice with the experimental lattice constant for both electrodes and exploit a recursive relation 
for the surface Green's functions of the form — a — jSgjS^ [l38|, where a and jS denote 



the diagonal and off-diagonal blocks of ES — H in the contacts, respectively. This equation is 
solved iteratively until g is converged, which corresponds to the materials-specific surface Green's 



4lh . we compute the surface GFs for 



function of a semi-infinite system. As in previous studies 11231 
all sites in a surface layer of two atomic planes. 

It is well known that the conductance is very sensitive to interference effects that arise form 
small atomic displacementllTi I29I] . These effects lead to weak oscillations in the total transmission 
at r = which are averaged over in most quantum transport measurements at room temperature. 
To account for this phenomenon we average t{E) over a small interval [Ef — A, Ef + A] around the 
Fermi energy Ef, whereas A = 50 meV ^ IksT to obtain a representative value of the zero-bias 
conductance for comparison with experiment. 

C. Extended Huckel-Hamiltonian 

In this investigation we have employed extended Hiickel theory (EHT), which has been widely 



58i 



59D, for the electronic structure 



used in molecular electronics calculations in the past 
calculations, but other semi-empirical methods can be used as well. Excellent semi-empirical pa- 
rameterizations exists for a wide array of systems, including metals, biomolecules (e.g. DNA) or 
semiconductors, but the limitations of this approach (e.g. regarding transferability) are also well 



known. In our calculations we employ a minimal basis set 



of non-orthogonal Slater type 
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orbitals |/), including s-, p- and d-wave functions for each atom. In this representation the matrix 
diagonal elements ha of the Hamiltonian H are approximated by the corresponding orbital ion- 
ization potentials, while the off-diagonal elements /i/y, / ^ j are chosen as {ha + hjj) • Sij/2. Here 
Sij = denotes the overlap matrix elements of S. The suitable choice of the basis functions 
permits a fast evaluation of the integrals The material specific orbital ionization potentials 
are extracted from experimental data or ab initio calculations. Here we employ a well established 
parameter set taken from Ref. [l62n . This choice of the Hamiltonian parameters is derived form 
electronic properties, and so it is partially suited for the description of the coherent transport prop- 
erties of a quantum system. In explorative calculations we found other semi-empirical parameter 
sets, which have been derived from structural properties and heats of formation (such as PM6), 
agreed less with prior theory and experiment. 
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FIG. 1: (Color online) Schematic representation of the different regions of a single molecule junction in 
the Landauer approach, (a) Definition of the extended molecule including a fraction of the electrodes. This 
permits the natural charge transfer 5n and screening effects as an organic system is attached to a metal 
surface, (b) Realization of a molecular junction by a phenyl-ring-based wire coupled via sulphur atoms 
to the [lll]-layers of gold leads. The semi-infinite leads are represented by their surface Green's function 
defined on the atoms of the dark shaded area. Horizontal lines indicate a possible division into a left and 
and a right part of the system, (c) Representative division of the system into principal layers to illustrate the 
recursive Green's function approach. 
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FIG. 2: (Color online) Convergence test of the conductance depending on the „principal layer"-thickness for 
several nano-junctions. (a) Model nanojunction of 45.2 A length and minimal cross section of one atom al- 
lowing for several different principal layer divisions, indicated by the marked lines below the conformation. 
On the right hand side the thickness of the principal layers of the actual division is indicated, respectively, w 
is given in units of the [111] atomic layer distance d\^\\\] = 2.35A (b) Corresponding conductance values for 
the upper described sets of „principal layers" for a silver and a gold contact, respectively. The dependence 
of the conductance on the layer division is also shown for similar metallic junctions with a minimal cross 
section of 2, 3 and 4 atoms, respectively, (c) Metal quantum wires with one conductance quantum, but 
with increasing length between 20 and 32 atomic layers show the same rapid convergence behavior with 
increasing principal layer thickness. 
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FIG. 3: (Color online) The total transmission t{E) of the model silver junction shown in Fig la, in the 
different levels of approximation, i.e. a changing number of „principal layers". Convergence with respect 
to the principal layer thickness is achieved as soon as w becomes is larger than 2. 
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FIG. 4: (Color online) The tilt- and twinning-angle dependence of the conductance of a silver point contact, 
(a) Only a moderate change of (G) is observed during tilting the electrodes up to 70 degrees, (b) The 
twinning of the electrodes between and 60 degrees results in a nearly constant conductance. 
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FIG. 5: (Color online) Representative examples of the generated conformations of silver nanojunction with 
an increasing number of surface vacancies. The presence of defects leads to a decrease of the conductance 
by up to 30%, indicated by the corresponding conductance values below the geometries. The number of 
vacancies in conformation is given in brackets. 
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FIG. 6: (Color online) Conductance values averaged over conformations with equal number of surface 
impurities. The error-bars indicates that the change in the conductance depends less on the number of 
defects, but more on their positions. 
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FIG. 7: (Color online) Transmission of silver nano-clusters. (a) Shows the top view on the examined nano- 
cluster conformations with 5, 7, 180, 220, and 260 atoms, (b) Cluster conformations energetically optimized 
on a silver substrate layer with a pyramidal electrode on top. (c) Calculated transmission function of the 
junction conformations shown in "b". The vertical line indicates the Fermi energy. 
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FIG. 8: (Color online) Configurations of the organic molecular wires studied in this work, (a) Oligopheny- 
lene molecules covalently bond to Aui9 clusters along the cristallographic [111] axis. One phenylene 
ring unit represents one principal layer, (b) Total transmission as a function of the energy of the shown 
oligophenylens with (dashed line) and without (sa 



lid line) the principal layer approximation in good qual- 
itative agreement with the DFT results of Ref. Il23n . In the layer approxiation one principal layer contains 
a single phenyle ring unit. The vertical line indicates the Fermi energy. {Inset) Length dependence of the 
conductance of the oligophenylene wires. The conductance decreases exponentially with the number of the 
phenyle rings in the wire in good agreement with experimental data. 



24 




time (ps) time (ps) 

FIG. 9: (Color online) Thermal influence on the conductance of a molecular wire at 300 K.(a) Fluctuation of 
the torsion angles 0i, 02 and 03 between the ring-units occuring in Au-h-R4, respectivlely. (b) Distribution 
histogram of the frequency of occurrence of a particular torsion angle, (c) Corresponding conductance 
(black) and average torsion angle (gray/turquoise) at the fluctuation process during 10 ps simulation 
time, (d) Zoom into the 4ps-5ps range, which shows that a short-time increase of the torsion angles (thin 
gray/colored curves, left axis) leads to a strong decay of the total transmission (bold black curve, right axis) 
of the nano wire. 

Tables 



TABLE I: Conductance of the molecular wires Rl, R2, R3, and R4 without and with the principal layer 
approximation compared to experimental results from Ref. \\52\ 
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molecule 


Gno layers (10 ^Gq) 


Gwith layers (10 ^Gq) 


Gexp (lO-^Go) 


Rl 


16.10 


14.60 


6.40 


R2 


4.10 


3.80 


1.16 


R3 


0.85 


0.80 


0.18 


R4 


0.20 


0.20 
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